Theory of superconductivity of carbon nanotubes and graphene 
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We present a new mechanism of carbon nanotube superconductivity that originates from edge 
states which are specific to graphene. Using on-site and boundary deformation potentials which do 
not cause bulk superconductivity, we obtain an appreciable transition temperature for the edge state. 
As a consequence, a metallic zigzag carbon nanotube having open boundaries can be regarded as a 
natural superconductor/normal metal/superconductor junction system, in which superconducting 
states are developed locally at both ends of the nanotube and a normal metal exists in the middle. 
In this case, a signal of the edge state superconductivity appears as the Josephson current which 
is sensitive to the length of a nanotube and the position of the Fermi energy. Such a dependence 
distinguishs edge state superconductivity from bulk superconductivity. 

PACS numbers: 74.20.Mn, 74.10.+V, 74.78. Na, 74.70.Wz 



Superconductivity in carbon nanotubes (NTs) has 
been attracting much attention due to its high supercon- 
ducting transition temperature, ^ 10 K. [l], 0] How- 
ever, it is well-known that superconductivity in low di- 
mensional (quasi- ID) systems is difficult to produce due 
to low density of states (DOS) 3], strong quantum fluc- 
tuations and other phenomena in such systems. More- 
over, metallic NTs exhibit ballistic transport properties 
at low temperatures, Q which suggests a weak electron- 
phonon (el-ph) interaction for the conducting electrons. 
It is surprising that superconductivity is realized in NTs 
at such high values of T^.. The mechanism of NT super- 
conductivity is a critical issue and determining it will be 
a valuable contribution not only to NT science but also 
to nanotechnology. 

Superconductivity has been observed in different types 
of NTs. Tang et al. reported Tc ~ 15 K for single-wall 
NTs (SWNTs) having a diameter of 0.4 nm. Q Takesue 
et al. found an abrupt drop in the zero-bias resistance 
at 12 K for multi-wall NTs (MWNTs) having an outer 
diameter of ~ 10 nm. [2] It is not straightforward to ex- 
plain the results obtained in experiments. For instance, 
the DOS at the Fermi energy of a SWNT appears to be 
too small to give rise to such high Tc- Kamide et al. § 
and Barnett et al. Q considered that curvature of (5, 0) 
SWNTs may increase the DOS. However, a large DOS 
may induce charge density wave (CDW) before super- 
conductivity occurs. Connetable et al. showed that (5, 0) 
and (3, 3) SWNTs undergo a CDW transition at temper- 
atures above room temperature. [2] Thus, small diame- 
ter SWNTs may be insulators. Moreover the curvature 
effect is negligible for MWNTs. The origin of NTs super- 
conductivity can not be explained by curvature-induced 
DOS and a new explanation is needed. 

Here, we focus our attention on the large local DOS 
(LDOS) given by edge states which are intrinsic to 



graphene. The edge states are electronic localized states 
that exist around the zigzag edge of graphene and a 
SWNT. [i, i] The energy dispersion of edge states is 
located near the Fermi energy [E-p = 0). The value of 
LDOS depends on the energy bandwidth (14^) of the edge 
states. Recent experiments involving scanning tunneling 
microscopy /spectroscopy (STM/STS) at the zigzag edge 
of graphene, [lOl [lH and angle-resolved photo-emission 
spectroscopy (ARPES) of Kish graphite, [l^ showed that 
the edge states are located below the Fermi energy and 
have a finite W . Although edge states of NTs have not 
been observed so far, it is possible to consider the edge 
states of a SWNT with a zigzag edge as well as those of 
a graphene sheet. 

In this letter, we calculate Tc as a function of W and 
E-p, and obtain an appreciable values for of the edge 
states of zigzag SWNTs and graphene. As a result, we 
predict that the superconductivity of a SWNT is given by 
a natural superconductor/normal metal/superconductor 
junction (SNS) system, in which superconducting states 
develop locally at both ends of the SWNT and a normal, 
ballistic state exists in the middle of the SWNT. Remark- 
ably, the bulk part of a SWNT need not be superconduct- 
ing since Josephson supercurrent flows in the middle as a 
result of the proximity effect when the superconducting 
edge states have different phases at both ends. We note 
that proximity-induced supercurrents have been observed 



in Ta/SWNTs/Au [13] and Nb/MWNTs/Al systems. [14| 
The Josephson current of a metallic zigzag SWNT de- 
pends on the length (L) and temperature (T). The 
amplitude of the current is proportional to exp(— L/^jv) 
when T < To where Tq ~ 20um/L K and ~ lO^K/T 
nm is the coherence length, [15| which is a characteris- 
tic feature of conventional SNS transport theory for the 
clean limit. [l3| A length dependence of the current dis- 
tinguishes edge-state superconductivity from bulk super- 
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conductivity. 

The edge-state superconductivity has the following 
advantages in explaining the experiment performed by 
Takesue et ai, Q (1) the edge states are robust against 
static surface deformation which is relevant for CDW in- 
stability, 17, [ll] (2) the el-ph interaction for the edge 



states is strong compared with that for delocalized states, 
and (3) Tc is sensitive to W and the energy position of 
Ep, which are all consistent with the fact that the su- 
perconductivity is sensitive to the junction structures of 
the Au electrode/MWNTs. [2] Enhancement of at the 
edge is important for understanding superconductivity 
of a general surface state, not only for graphite materials 
but also for noble-metals such as gold. jl9i] 

The edge states are zero-energy {E{k) = 0) eigenstates 
of the nearest-neighbor (nn) tight-binding Hamiltonian, 
Hnnl'^ik)) = E{k)\'^{k)). 0, [II] fc = k • ai is the 
wavevector around the tube axis where ai is the unit 
vector along the edge (see Fig. [UJa)). The edge states 
exist for 2tt/3 < k < 47r/3. The wavefunction is written 
as 



(1) 



where |0(Ri)) is the 2pz orbital, and Ci{k) is the ampli- 
tude at Hi and has a value on one of the two sublattices 
(A and B) of graphite. Q In the direction of the SWNT 
axis, the magnitude of Ci (k) quickly decays from the edge 
to the interior region. The localization length is given by 
^{k) = -|T|/21n|2cos(fc/2)| (27r/3 < A: < 47r/3) where 
T = 2a2 — ai is the translation vector, [l^, [2l| When 
we incorporate the next nearest-neighbor (nnn) trans- 
fer integral 7^ into the Hamiltonian, the energy disper- 
sion of the edge states becomes E{k) — 7„(2cosfc + 1) 
(27r/3 < k < An/S) where the value of 7„ = 0.3 eV 
is adopted. [20l . [2a | The calculated results explain the 
STS [lElili and ARPES [13] experiments. Hereafter we 
treat W (= 7n) and the position of the Fermi energy as 
independent parameters. 

The el-ph interaction for the edge states shows a dif- 
ferent behavior from that for delocalized states. The el- 
ph interaction consists of on-site and off-site deformation 
potentials. 23| It is pointed out that, for a backward 



scattering of delocalized states, the on-site deformation 
potentials on two sublattices cancel with each other due 
to a phase difference of the wavefunction at the two sub- 
lattices. 2^ This is a reason why metallic NTs show a 
ballistic transport property. However, the cancellation of 
the on-site deformation potential does not work for the 
edge states since the wavefunction of the edge state has 
an amplitude only on one of the two sublattices. Fur- 
thermore, because of a lack of translational symmetry at 
the edge, a strong el-ph interaction for optical phonon 
modes is expected for the edge. Thus the understanding 
of the el-ph interaction for the edge states is essential for 
the present problem. 



The el-ph Hamiltonian is defined by 7^ = Hq + 
Hint, where Ho represents the edge states and phonon, 
and Hint is the el-ph interaction. Hq is given by 

J2k^i^)4:'^k + Z]q,i.'^i'(q)^q,i.^q,i'; whcrc Cfc is the an- 
nihilation operator of edge state and bq^i, is the annihi- 
lation operator of i^-th phonon mode with momentum 
q and energy uj^{q). For a graphite unit cell, there are 
six phonon eigen- modes; out-of- plane tangential acous- 
tic/optical mode (oTA/oTO), in-plane tangential acous- 
tic/optical mode (iTA/iTO), and longitudinal acous- 
tic/optical mode (LA/LO). uj^{q) and phonon eigenvec- 
tor are obtained by solving a 6 x 6 dynamical matrix. [2^ 
The el-ph interaction is given by 



Hint — 



k,k' qt,v 



where Nu_ is the number of graphite unit cells in a SWNT, 
and a^j,, (q) is the el-ph coupling connecting two edge 
states k and k' by z/-th phonon mode with momentum 
q. Due to the momentum conservation along the edge, 
k' — k+q (q = q-ai), while the wavevector perpendicular 
to the edge qt (= q-T) is needed to sum over the Brillouin 
zone. 

We calculate a'^/^, (q) using the deformation potential, 
SV — —J2p'^^i'^p) ' u(Rp), where u(Rp) is the dis- 
placement vector and w(Rp) is the pseudo-potential of 
a carbon atom at Rp. The present pseudo-potential 
is used for calculating resonance Raman intensity in 
which the calculated results explain chirality and diame- 
ter dependence of Raman intensity quantitatively. (23| 
u(Rp) can be expanded by phonon normal modes as 

u(Rp) = Eq,.(^''(q)/V2^r)(feq,.+fe-q,.)e-(Rp)e''^•I^^ 

where eq(Rp) is the normalized eigenvector at Rp 



and A'^{q) — fij \J •mc<^i,{p^ is the phonon ampli- 
tude. From (^'(fc')|(5y|*(fc)), we obtain ^^/(q) = 
A''(q)M^j.,(q)/V2, where M^j.,(q) is the el-ph matrix 
element defined by 

M,-,,(q) ^ -E(*(fc')|VKRp)|*(fc)) •e:;(Rp)e'^'i-^^. 

p>0 

(3) 

Putting Eq. ^ to Eq. we see that M4,(q) 

consists of the on-site (0(Ri)| Vu(Rp)|0(Ri)) and off- 
site ((/)(Rp)|Vu(Rp)|(/>(R,;)) atomic deformation poten- 
tials (Rp 7^ Rj). The off-site atomic deforma- 
tion potential does not contribute to M^^,, (q) because 
((/)(RjJiyi;(Rp)|0(R,)) vanishes for |RpeA - R*| > 
|ai|. [iil We note that Eq. ^ includes the effect of 
boundary. To show this, we illustrate several carbon 
atoms (Rp>o) near the zigzag edge and a fictitious atom 
(Rp<o) in Fig. [T{a). The on-site deformation poten- 
tial at Ri is given mainly by the vibrations of carbon 
atoms at R2 and R3 as — X]p=2 3(</'(R-i)l^^(r^p)l</'(I^i))' 
eq(Rp)e''''^p. This on-site deformation potential would 
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be canceled by the fictitious carbon atom at R_i since 
eq(Rp) {p = —1,2,3) points at the same direction. 
Namely, the on-site deformation potential is enhanced 
at the edge. This enhancement may be a reason why the 
tunnel current is unstable at the edge, [lot 

For a (n, 0) SWNT, k for the edge states becomes dis- 
crete as k{i) = 2tt/3 + 2Tri/n (i = 1, . . . , n/3 — 1) due to 
the periodic boundary condition around the axis. We de- 
note an edge state by the integer i and write .^^^^.^ (q) 
as Mj^(q) for simplicity. Putting k{i) to ^{k), we ob- 
tain ^(fc) < dt/2 where dt = n|ai|/7r is diameter for the 
SWNT. 

In Fig.IHb) and (c), we plot \M[i{q) \ and |M|'g(q)|, re- 
spectively, for the (60,0) SWNT {dt « 5 nm). \M^i{q)\ 
is chosen as an example that ^(1) ^ 22A and ^(4) ~ 
4.4A are much longer than the carbon-carbon bond 



T,{k,iujn) is the self-energy, which is determined self- 
consistently by 



while |Mgg(q)| is chosen as another 



length Occ 1.4 A, 

example that ^(6) ~ 2.4A and ^(9) ^ 0.9A are compara- 
ble to flee- As for acoustic modes, the LA mode couples 
strongly to the edge state. The oTA mode contributes 
to |M{'4(q)|, whereas the iTA mode is negligible. Since 
the LA and oTA modes change the area of a hexagonal 
lattice, they contribute to on-site deformation potential. 
For optical modes, the iTO and LO modes are important. 
|M{4°(q)| decreases with increasing qt, while jMjlP (q)| 
increases with increasing qt. As shown in Fig. W^c), the 
deformation potential is stronger for the smaller local- 
ization length. The behavior of the iTO and LO modes 
is due to the boundary deformation potential. To prove 
this, we show in the inset, the matrix element without 
the boundary, which is defined by Eq. ([3]) including (ficti- 
tious) carbon atoms at p < in Fig. [11(a) . The boundary 
deformation potential depends on the direction of and 
the magnitude is maximum when ej^ is parallel to T. ej^^*-* 
(eq°) has a large element parallel to T when qt < \/3q 
{qt > V3g) as shown as a vertical line in Fig.[T](b) and (c). 
On the other hand, e°"'"'^ is perpendicular to the SWNT 
axis (or parallel to ai x T) and the boundary effect of 
the oTO mode does not appear in |M°^°(q)|. 

Now we apply |a^j,/(q)| to the Eliashberg equation. 
The Eliashberg equation includes the effects of phonon 
retardation and electron self-energy^ which are not taken 



into account in the BCS theory. [25|, |26[ The phonon re- 
tardation is included by the Matsubara frequency: Un = 
kBT{2n + l)7r where n is integer and \ujn\ < ujd where 
tJD = 0.2 eV is the Debye energy. Q Since the gap func- 
tion, A{k,iujn), vanishes at Tc, the Eliashberg equation 
can be linearized at Tc to get the gap equation: 



A(fc, iujn) 



2kBTc 



^ |<fc,(q)pc..(q) 

Nu (a'„-tJ™)2+a)3(q) 

X |G(fc',iw™)|2A(fc',iw™), 



(4) 



S(fc, iw„) 



2kBT, 



<,,(q)pc.,(q) 

^ - UJmY +^1{Ol} 



G{k',iuJrn)- 

(5) 



where G{k,iujn) is a thermal Green function of electron, 
G{k,iujn) = {it^n — {E{k) — Ep) — Y,{k,iuj„))~^ . Here, 



After calculating 'S,{k,iujn) in Eq. ([5]), we solve Eq. (jl]). 

In Fig. [21 we show Tc as a function of W for (n, 0) 
SWNTs with n — 30, 60, and 90, where we assume = 
—W 12. Tc decreases with increasing W and Tc vanishes 
at critical values, W^. The increase of W corresponds to 
the decrease of the LDOS around the Fermi energy. The 
values of Wc are 0.46 eV and 0.37 eV, respectively for n = 
30 and n > 60. Those values of Wc are close to 7„ (the 
dashed line in Fig.[21). When n {dt ) is relatively small, all 
edge states couple strongly to the boundary deformation 
potential since max(^) ^ dt/2. The strong el-ph coupling 
for n = 30 makes Wc larger than that for n > 60. It is 
also noted that excluding optical modes makes Tc and Wc 
both smaller. In this case, we obtain Tc ~ 70 K at W = 
eV and Wc 0.21 eV for n = 60. Although the values of 
A'^(q) for optical modes are smaller than those of acoustic 
modes, the iTO and LO modes contribute to Eqs. (|4|) 
and ([U because of the large values of |Mj^ (q)| due to 
the boundary deformation potential. A large value of n 
corresponds to the zigzag edge of graphene. Remarkably, 
Tc for n = 120 has a curve quite similar to Tc for n — 
90. This suggests that Tc converges and n = 90 is large 
enough to represent a graphene. 

It is important to note that the calculated Tc is sen- 
sitive to the energy position of Ep. We plot Tc as a 
function of Ep for a (60,0) SWNT with W = 7„ in the 
inset of Fig. [2] When Ep exists at the top of the energy 
band of the edge states, Tc becomes less than 1 K. When 
Ep ^ —0.1 eV, Tc decreases rapidly since the inelastic 
scattering process is suppressed by the absence of the 
scattered state. This is a reason why Tc is sensitive to 
the position of Ep. We also calculated Tc for extended 
states around the Fermi energy of the (60, 0) SWNT us- 
ing the Eliashberg equation. The calculated Tc is less 
than 0.1 K. Thus the extended states do not contribute 
toTc. 

The observed Tc should be smaller than our estimation. 
In fact, a lattice defect along the edge decreases LDOS 
and reduces Tc. The Coulomb repulsive interaction might 
decrease Tc, too. Fujita showed that the edge states 
develop a local ferro-magnetism in the presence of a large 
Hubbard U comparable to W. Q Since the edge states 
are localized at the edge, they might have a quantum fluc- 
tuation intrinsic to ID system. The Tomonaga-Luttinger 
liquid theory may be suitable to calculate the correlation 
function. 

In summary, using the Eliashberg equation, we clarify 
that W and Ep position is sensitive to Tc of the edge 
states in SWNTs and graphene. The rather high value 
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FIG. 1: (Color online) (a) Carbon atoms at Rp>o near the zigzag edge, and a fictitious carbon atom at R_i are illustrated to 
show an enhancement of the on-site deformation potential at the boundary. The boundary deformation potential is large for 
optical modes whose is parallel to T. (b,c) \M^j{ci)\ of the (60,0) SWNT with (b) fc(l)(= 77r/10) k{4){= 4tv/5) and (c) 
fc(6)(= 137r/15) — > fc(9)(= 297r/30) are plotted as a function of qt where q — tt/IO. The inset is the matrix element including 
(fictitious) carbon atoms at Rp<o. Three solid/dashed curves represent acoustic/optical phonon modes: oTA/oTO (green), 
iTA/iTO(blue) and LA/LO (red). The vertical dashed lines represent qt = VSq {VS = lTl/|ai|). 
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FIG. 2: The dependence of Tc on W is plotted for (30,0), 
(60, 0) and (90, 0) zigzag SWNTs. The Fermi energy is as- 
sumed to be located at the center of the band {Ef = — W/2). 
The nun hopping gives W = 0.3 eV. (inset) The dependence 
of Tc on -Ef for the (60, 0) SWNT with W = j^. 



of Tc obtained is a result of LDOS enhancement by the 
edge states, and the on-site and boundary deformation 
potentials of the el-ph interaction for the edge states. If 
nanotube superconductivity is given by el-ph interaction, 
the edge-state superconductivity is a unique candidate 
since Tc of the bulk is negligible. Edge (surface) state su- 
perconductivity is potentially a key concept for designing 
superconductors on the nanometer scale. 
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